Jump to content

User:Compsonheir/Stefan problem

From Wikipedia, the free encyclopedia

In applied mathematics, the Stefan problem constitutes the determination of the temperature distribution in a medium consisting of more than one phase, for example ice and liquid water. While in principle the heat equation suffices to determine the temperature within each phase, one must also ascertain the location of the ice-liquid interface. Note that this evolving boundary is an unknown (hyper-)surface: hence, the Stefan problem is a free boundary problem.

The problem is named after Jožef Stefan, the Slovene physicist who introduced the general class of such problems around 1890, in relation to problems of ice formation. This question had been considered earlier, in 1831, by Lamé and Clapeyron.

Mathematical formulation

[edit]

Consider a substance consisting of two phases which have densities ρ1, ρ2, heat capacities c1, c2 and thermal conductivities k1, k2. Let T be the temperature in the medium, and suppose that phase 1 exists for T < Tm while phase 2 is present for T > Tm, where Tm is the melting temperature. With a source Q, the temperature T satisfies the partial differential equation

in the region where T < Tm,

in the region where T > Tm. However, there is an additional equation determining the location of the interface between the two phases. Let V be the velocity of the phase boundary; we adopt the convention that V is positive when the phase boundary is moving towards phase 2 and negative when moving towards phase 1. Then V is determined by the equation

.

Here the subscript 1 is a shorthand for

and vice versa for 2.

Derivation

[edit]

The Stefan condition

is the key component of this problem. It can be derived in a similar fashion to the Rankine-Hugoniot conditions for conservation laws, which we do below.

First, we note that the internal energy of the material is

where L is the latent heat of melting. The differential equations and Stefan condition are consequences of the fundamental relation

where k is the thermal conductivity and Q represents diabatic heating.

Unfortunately, the differential equation for the internal energy does not make sense in a neighborhood of the region where the substance is melting: the internal energy is discontinuous due to the additional latent heat necessary to melt a solid, and the thermal conductivity is discontinuous too.

To remedy this situation, we interpret the equation for the internal energy as a conservation law by integrating the PDE over a small box in space-time containing a point on the surface

using the divergence theorem and then taking the limit as the size of this box shrinks to zero. For the sake of an unambiguous definition we take

The vector n points in the same direction whether the limit is taken from the solid or liquid side of the phase boundary, since the gradient always points in the direction of increasing T. Now, note that the jump of E going from the solid to liquid phase is

See also

[edit]

References

[edit]
[edit]

Category:Partial differential equations